Tunneling spectra of layered strongly correlated d-wave superconductors 
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Tunneling conductance experiments on cuprate superconductors exhibit a large diversity of spectra 
that appear in different nano-sized regions of inhomogeneous samples. In this letter, we use a mean- 
field approach to the tt't" J model in order to address the features in these spectra that deviate 
from the BCS paradigm, namely, the bias sign asymmetry at high bias, the generic lack of evidence 
for the Van Hove singularity, and the absence of coherence peaks at low dopings. We conclude that 
these features can be reproduced in homogeneous layered d-wave superconductors solely due to a 
proximate Mott insulating transition. We also establish the connection between the above tunneling 
spectral features and the strong renormalization of the electron dispersion around (0, n) and (n, 0) 
and the momentum space anisotropy of electronic states observed in ARPES experiments. 
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Introduction - Tunneling and angle-resolved photoe- 
mission spectroscopy (ARPES) are unique experimen- 
tal techniques in that they probe the single electron 
microscopic physics of strongly correlated systems like 
the cuprate <i-wave superconductors (dSC). Even though 
both techniques support the presence of well defined Bo- 
goliubov nodal quasiparticles in these materials, away 
from the nodes deviations from the BCS paradigm are 
encountered P,SIESSIE0llil3- Simultaneously de- 
scribing these deviations in ARPES and tunneling exper- 
iments is of utmost importance to understand the nature 
of the underlying strong local correlations. 

In this letter, we focus on the differential conductance 
measured by scanning tunneling microscopy (STM) ex- 
periments on the cuprates. In particular, we address 
three key aspects observed in the data as the electron 
concentration is increased toward half-filling, namely, (i) 
an increasing asymmetry of spectra for large positive and 
negative bias; (ii) the ubiquitous absence of peaks caused 
by the underlying Van Hove (VH) singularities in the 
quasiparticle density of states (DOS); (Hi) the gradual 
depletion of the superconducting (SC) coherence peaks. 

First, we argue that BCS expectations for dSC are 
at odds with the above features. BCS theory generally 
predicts that tunneling spectra of quasi-2D SC materi- 
als display VH singularity peaks in addition to two SC 
coherence peaks 0,0. However, these VH peaks are 
not observed by expe riments in different cuprate families 
[Bi2212 HE El El El, B12201 Q, LSCO El, YBCO 
Notably, even the SC coherence peaks are absent 
in certain tunneling spectra U& u& Uj| • 

We then show that the aforementioned non-trivial 
spectral features are reproduced in homogeneous SC sys- 
tems close to the Mott insulating state. Specifically, 
we use a mean-field (MF) theory of doped Mott insu- 
lators [lil Ej| that accounts for many fingerprints of 
strong correlation physics in ARPES data, including the 
strong band renormalization close to (0, tt) and (it, 0) [13] 
and the sharp momentum space anisotropy of electronic 



states, known as the nodal-antinodal dichotomy 
which we show are intimately connected to STM data. 
We thus obtain a consistent theoretical microscopic de- 
scription of ARPES and tunneling spectroscopy. More- 
over, we explicitly show that the above deviations from 
BCS spectra follow from how the local electron Coulomb 
repulsion and the resulting short-range antiferromagnetic 
(AF) correlations affect the phenomenology of dSC. 

Tunneling spectra of BCS d-wave superconductors - 
We consider the tunneling across a normal-metal - in- 
sulator - superconductor junction in the case of a 2D 
superconductor. Assuming specular transmission across 
a thin planar junction within the elastic channel, the dif- 
ferential tunneling conductance perpendicular to the 2D 
layers can be essentially equated to the quasiparticle DOS 
|2ll \22\ and, at zero temperature, it is given by |23| 



dl 



J- = 4eirM 2 N(E F ) ^ A(V, k) 



(1) 



where M is the tunneling matrix element, N(Ep) is the 
normal-metal DOS, V is the bias of the SC sample with 
respect to the metal and A(u>, k) is the electron spectral 
function in the SC sample. The sum in Q is over the 
2D layer's Brillouin zone. 

In Fig. ^ we use O together with the well known re- 
sults for the BCS electron spectral function to plot 
the BCS tunneling differential conductance. In particu- 
lar, we take the normal state dispersion s k determined 
by the 1 st , 2 nd and 3 rd nearest- neighbor (NN) hopping 
parameters t, t' and t" with t' = -2t" = -0.12t [2J so 
that the normal state Fermi surface lies close to (0, 7t) 
when the electron density is 1 — x — 0.80 (inset of Fig. 
EJ, as observed in LSCO [jfj Bi2201 [jj and in the an- 
tibonding sheet of Bi2212 |28j. In addition, we consider 
the d-wave electron pairing gap function = A cos(26 l ) 
with A = O.li. The resulting dl/dV spectrum [see dash- 
dot line in Fig. da)] depicts two sharp coherence peaks 
at V = ±A due to the SC quasiparticle energy disper- 
sion E k = [el + A 2 ,] 1 / 2 saddle-point at the SC gap en- 
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FIG. 1: BCS tunneling differential conductance for t' — 
-2t" = -0.12/; and: (a) A = O.lt with x = 0.20 (dash-dot 
line) and x = 0.10 (solid line); (b) A = 0.2/j, x = 0.10 with 
(solid line) and without (dash-dot line) Lorentzian broadening 
F = O.Olt. Bias is in units of J = t/3. Plots are obtained for a 
2000 x 2000 mesh of the Brillouin zone and energy resolution 
of0.0067t Inset: Fermi surface of e k with t' = -It" = -0.12t 
for x = 0.20 (dash-dot line) and x — 0.10 (solid line). 



ergy. This behavior is consistent with experimental spec- 
tra around optimal doping H S H [H El . 

However, if the electron concentration is increased to 
1— x = 0.90 the normal state Fermi surface intersects the 
(0, 7r) — (tt, 7t) line at a point, which we denote by kA, that 
moves away from (0,7r) (inset of Fig. In this case, 
the VH singularity in Ek at (0, tt) and (n, 0) shows up as 
a separate peak in the BCS tunneling spectrum at V — 



-E- 



VH 



-fe=(0,7r) 



+ A 2ji/2 [pig. I^a)]. This additional 
peak appears behind the negative bias SC coherence peak 
unless k,A ~ (0, tt) or Eyn — A is so small that spectral 
broadening smears the VH and the SC coherence peaks 
into a single feature. However, as shown in Fig. ^b), the 
above double peak structure still survives if, as expected 
for real materials upon increasing electron concentration, 
the d-wave gap is doubled to A = 0.2i and a broadening 
width T = 0.01* is included H|. 

This behavior contrasts with the experimental obser- 
vation that no such double peak structure occurs as the 
hole doping level is decreased [1 H 0, 0, Q and im- 
plies that the actual dispersion along (0, 7r) — (it, tt) in 
the cuprate materials is of the order of the SC gap en- 
ergy scale, namely ~ O.li, which is an order of magni- 
tude lower than t. Such a flat dispersion in the antinodal 
region is supported by ARPES data in Bi2201, Bi2212 
and YBCO [2(j. Since ARPES observations also con- 
strain the much larger nodal dispersion energy scale, as 
well as the curvature of the Fermi surface, several fine- 
tuned phenomenological parameters are required to fit 
the observed normal state dispersions [Fig. |2fa) displays 
the dispersions obtained by Norman et al. from experi- 
mental fits involving up to five hopping terms |30l l3l|]. 
As Ref. points out, the ubiquitous discrepancy be- 
tween the energy scale of the nodal dispersion and that 
of the extended flat region around (0, tt) throughout the 
various cuprate families and doping levels suggests that 



interactions strongly renormalize the electron dispersion. 
Below, we argue that such phenomenological dispersions 
are natural to doped Mott insulators. 

Renormalization of single electron properties - To dis- 
cuss how the single electron dispersion relation and 
the momentum space distribution of spectral weight are 
renormalizcd in doped Mott insulators we consider the 
tt't"J model, whose numerical calculations reproduce 
many spectral features of the cuprates [H 0| . 

In the half-filled AF insulator, the electron dispersion 
is renormalized by the strong local AF correlations. In 
particular, NN hopping is heavily frustrated by the AF 
background and, for t' = t" = 0, holes hop coherently 
within the same sublattice by virtue of spin fluctuations 
[35| . As a result, the width of the nodal dispersion is 
reduced by a factor of ~ 10 from the bare value 8i 
down to 2. 2 J 33] . t 1 « — 2t" control the dispersion 
along (0, it) — (tt, 0) and, since they describe intrasub- 
lattice hopping processes which are not frustrated by AF 
correlations, the dispersion width along this line is only 
renormalized by a factor ^2 — 3 [33| . 

Away from half-filling, holes distort the AF spin back- 
ground and change the surrounding spin configuration 
in order to coherently hop between NN sites and thus 
gain extra kinetic energy. Interestingly, exact diagonal- 
ization studies show that the resulting local spin corre- 
lations strongly renormalize t' and t" j3||. This effect, 
which is predicted in spin liquids |3?|. is consistent with 
the experimentally observed doping induced flatness of 
the dispersion in the antinodal region |2f3 | since the dis- 
persion along (0,7r) — (7r,0) is controlled by t' and t". 

The "doped carrier" MF approach introduced in Refs. 
accounts for the above interplay between the ex- 
change energy of localized spins and the kinetic energy 
of delocalized holes. Here, we present only a qualitative 
description of the MF theory and defer the readers to the 
above references for formal details. This theory describes 
doped Mott insulators in terms of two different fermions: 
(i) "dopons" , which have the same charge +e and spin- 
1/2 as holes, describe vacancies surrounded by an AF 
spin configuration; (ii) "spinons" , which have no electric 
charge, describe spin-1/2 excitations of the spin back- 
ground. The MF Hamiltonian thus has four fermionic 
bands: e^ fe = -e+ fc and q k = -e+ k [Fig. I^b)]. If 
dopons and spinons do not hybridize the high energy 
bands e^ k describe the dynamics of holes in an AF back- 
ground and, at MF level, the spinon ef k bands do not 
contribute to electron spectral properties. The process 
of hybridizing dopons and spinons, which describes the 
change in spin correlations due to hole hopping, couples 
spin and charge dynamics and transfers electron spec- 
tral weight from e^ k to the lower energy ef k bands [Fig. 
EJc)]. This spectral weight transfer is not uniform in 
momentum space and reflects the dispersion of the do- 
pon band, which has lower energy at (tt/2,tt/2) than 
at (0,7r). As a result, the ei & band has more spectral 
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FIG. 2: Energy dispersions along (0, 0) — (tt, tt) — (0, tt) — (0, 0). 
(a) Normal state disp ersions using parameters from Ref. |30(| 
(solid line), Ref. 31. (dash-dot line) and e k defined in main 
text (dashed line), (b) MF dispersions ej fc and e^ fe for x = 
0.18 and parameters used in Ref. llSl (c) MF spectral function 
with Lorentzian broadening F — J/5 where J = 150meV. 
Dashed line in (b) and (c) plots [(e%) 2 + A 2 k ] 1/2 x sgn(e^). 
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FIG. 3: MF tunneling differential conductance for (a) x = 
0.20, (b) x = 0.10 and (c) x — 0.05. (d) Same spectra with 
broadening given by £"(cj) = —lo 2 /(5J). Momentum space 
and energy resolution are the same as in Fig. Q Insets of 
(a)-(c) plot minimal gap locus of tf k - Vertical bars in (a)-(c) 
denote VH singularity and SC gap energies. 



weight near the nodal points than in the antinodal region 
|l8| . in agreement with the nodal-antinodal dichotomy 
in ARPES data 0, H, u| I n conformity with other ap- 
proaches to the related Hubbard model [33] , we conclude 
that this dichotomy results from short-range correlations, 
namely AF correlations 0,H^|, in dSC. 

The low energy dispersions ef k derive from the spinon 
bands and, thus, from the spin exchange interaction. 
Therefore, they are controlled by a single energy scale, 
namely J, which leads to an almost flat dispersion near 
the antinodal points. We thus find that the spin back- 
ground strongly renormalizes the electron dispersion in 
the antinodal region. To explicitly compare the MF dis- 
persion with the experimental fits to ARPES data we 
introduce the underlying normal dispersion such that 
[(e^) 2 + A 2 .] 1 / 2 equals: (i) e+ fc at k A , (0, tt) and at the 
nodal point; (ii) e+ fc at (0,0) [Figs. Hb) and^c)]. In 
Fig. 01a) we plot ej^ for x = 0.18 against the normal 
dispersion fits from Refs. l3dL l3ll along major symmetry 
directions. Clearly, below the Fermi level e k captures 
the energy scale of the nodal dispersion and the flatness 
around antinodal points. ARPES does not probe the dis- 
persion above the Fermi energy and both experimental 
fits use band calculation results to fix the energy at (tt, tt) 
|30| . whence the mismatch between e k , obtained within 
the tt't" J model context, and these fits above the Fermi 
level. Such a mismatch is supported by exact diagonal- 
ization calculations of the tt't" J model which indicate 
that the dispersion above the Fermi level is less disper- 



sive than expected from bare hopping parameters |34j . 

Tunneling spectra of strongly correlated d-wave super- 
conductors - We now consider the tunneling conductance 
for the above MF theory of doped Mott insulators. We 
use the same parameters as in Ref. ^| to obtain the 
dl/dV curves at x = 0.20, x = 0.10 and x = 0.05 [Figs. 
|3fa)-|2fc)]. We point out that the MF theory produces a 
sharp electron spectral peak at the e^ fc band. However, 
we expect this spectral peak to be significantly broad- 
ened if one goes beyond the MF approximation since ef k 
are high energy bands which have many channels to de- 
cay. Hence, in order to better compare to experiments, 
in Fig. OHd) we depict the above spectra with a frequency 
dependent broadening given by 2"(a->) = — lo 2 /(5J). 

Looking at Figs. Et a )"EI c ) we recognize the separate 
contribution to the tunneling DOS from bands e^ fc , 
and e^ fe . The positive energy ej fc band does not ap- 
pear since it has vanishing spectral weight [see Fig. EIc)]. 
The presence of the negative energy bands e 2 k and k 
leads to a peak-dip- hump in the dl /dV curve as observed 
by experiments ja, Il2| . This negative bias hump, which 
moves to higher energy as doping is lowered 0, 0] , de- 
rives from the lower Hubbard band and renders the tun- 
neling conductance plots highly asymmetric for large pos- 
itive and negative bias. At smaller values of x the spec- 
tral weight in the e^ k and e^~ fc bands is transfered to the 
high-energy hump and, thus, the above bias sign asym- 
metry is enhanced upon underdoping as seen in STM 
data 0. [3. 0. 0. In particular, the ratio of integrated 
spectral weight at positive bias over that at negative bias 
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is 2x/(l—x), in perfect agreement with sum rules applica- 
ble to the generalized-^ J model |39j ■ The above behavior 
is a hallmark of models where strong local Coulomb re- 
pulsion between electrons makes it easier to add a hole 
to the sample than to add an electron. 

Figs. Et a )-Et c ) depict the doping evolution of the min- 
imum gap locus, which reproduces experiments [2fij . and 
denote the VH singularity and SC gap energies, Evh and 
A respectively, by the small vertical bars. Changing from 
x = 0.20 to x = 0.10 the wave vector fc^ distances from 
(0,7r), yet Evh — A remains extremely small since the 
dispersion is very flat in the antinodal region. Therefore, 
no extra peak develops next to the SC coherence peak, 
unlike the cases depicted in Fig. ^ The energy difference 
Evh — A increases for x = 0.05, yet, we observe a strong 
suppression of the SC coherence peak in the tunneling 
spectrum of Fig. |3fc) instead of a double peak struc- 
ture. Below, we explain this behavior, which is in stark 
contrast with weak coupling BCS predictions. 

The curves in Fig. [3] show how the SC coherence 
peak evolves with doping level. The x — 0.20 spec- 
trum displays well defined SC coherence peaks, in agree- 
ment with spectra of samples around optimal doping 
H H H [H [H 03. In the underdoped x = 0.10 spec- 
trum the SC coherence peaks are pushed to higher energy 
due to the SC gap doping dependence and lose spectral 
weight, as seen in experiments [U01H>E1- The spec- 
trum for the SC state at x = 0.05 hardly resolves the SC 
coherence peaks as found in deeply underdoped Bi2212 
0, and NaCCOC 0, E3- The gradual depletion of 
the SC coherence peak intensity as the half-filled state is 
approached follows from the correspondingly gradual de- 
pletion of spectral weight in the antinodal region which 
underlies the formation of the low energy spectral weight 
arcs observed both in ARPES experiments 0, IM an( i 
in theoretical approaches 0, |33- In our calculation, 
this behavior explicitly results from the presence of short- 
range AF correlations in dSC. 

As we emphasize in Fig. Old), the doping evolution 
of the dl/dV curves leads to qualitatively different line 
shapes around optimal doping and in the deeply under- 
doped regime. We refer to Fig. 2 of Ref. |a to show 
that the above results reproduce experimental evidence 
for different types of spectra in SC cuprate materials. 
Since this STM spectral diversity occurs in samples which 
are inhomogeneous on the nano-scale 0, 0, la, [13, |4fJ, |4JJ , 
it has been proposed that these non-trivial line shapes 
are intrinsically related to the presence of the underly- 
ing sample inhomogeneity |4(lj . The possible relevance 
of coexisting orders and phase separation scenar- 
ios 0, 0| has also been pointed out in the literature. 
However, here we show that the aforementioned unusual 
and diverse spectra can be reproduced in homogeneous 
SC systems solely due to the presence of strong correla- 
tions that follow from a proximate Mott insulating state. 
This theoretical assertion supports that electron dynam- 



ics in these materials is, to a large extent, determined 
by properties which are local even when compared to the 
nano-sized sample inhomogeneities, in conformity with 
experimental results [4l| . 
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